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ABSTRACT 



Aims. We carried out a radial-velocity survey to search for planets around metal-poor stars. In this paper we report the discovery of 
two planets around HIP 1 1952, a metal-poor star with [Fe/H]= -1.9 that belongs to our target sample. 

Methods. Radial velocity variations of HIP 1 1952 were monitored systematically with FEROS at the 2.2 m telescope located at the 
ESO La Silla observatory from August 2009 until January 201 1. We used a cross-correlation technique to measure the stellar radial 
velocities (RV). 

Results. We detected a long-period RV variation of 290 d and a short-period one of 6.95 d. The spectroscopic analysis of the stellar 
activity reveals a stellar rotation period of 4.8 d. The Hipparcos photometry data shows intra-day variabilities, which give evidence 
for stellar pulsations. Based on our analysis, the observed RV variations are most likely caused by the presence of unseen planetary 
companions. Assuming a primary mass of 0.83 M Q , we computed minimum planetary masses of 0.78 Mj up for the inner and 2.93 M Jllp 
for the outer planet. The semi-major axes are a\ = 0.07 AU and ai = 0.81 AU, respectively. 

Conclusions. HIP 11952 is one of very few stars with [Fe/H]< -1.0 which have planetary companions. This discovery is important 
to understand planet formation around metal-poor stars. 

Key words, star: general - star: individual: HIP 1 1952 - planetary system - technique: radial velocity 



1. Introduction 

Current results of the exoplanet surveys strongly suggest a cor- 
relation between a star's stellar metallicity and its probabil- 
ity of hosting planets, in particular for main-sequence stars 
(e.g., IFischer & Valentil 2005; I Johnson et al. 2010b . According 
to these studies, the detection rate of planets decreases with 
metallicity. However, the conclusions of IFischer & Valentil 
(2005) might be affected by an observational bias, since these 
authors did not have similar numbers of stars in their survey per 
bin of metallicity. Therefore, it is crucial to understand if either 
the high stellar metallicity triggers planet formation or the metal 
enhancement of stars is caused by the formation of planets. 

In the last years, exoplanet surveys tried to bridge this 
gap, starting to include more metal-poor stars in their samples. 
|Sozzetti et al. (2009)| conducted a three-year RV survey to look 
for planets around metal-poor stars down to [Fe/H]= -2.0 and 
found no evidence for short-period giant planets within 2 AU 
from the central star. Santo s~etal. (201 1)| performed a similar 
survey, only down to [Fe/H] — - 1 .4 and found three moderately 



metal-poor stars hosting a long period giant planets (P> 700 d). 
A hot Saturn and a hot Jupiter have been found transiting two 
moderate metal-poor stars, with [Fe/H]= -0.46 and -0.4, re- 
spectively (Bouchy et al. 2010 Simpson et al. 201 1 1. 
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In June 2009 we started a survey to search for planets 
around metal-poor stars. The target sample includes 96 metal- 
poor A and F stars. Our target list includes stars with metallic - 
ities in the range -4.0 <[Fe/H]< -0.5. As part of this work, 
|Setiawan et al. (2010)| found a planet around an extremely metal- 
poor red horizontal branch star with a short period of 16.2 d. 
We notice that its [Fe/H]= -2.1 is not included in the metal- 
licity range covered by the surveys of Soz zetti et al. (2009)| and 
ISantos et al. (201 1)| 

These recent observations have started to disclose the realm 
of planets at rather low stellar metallicities, indicating that metal- 
licity may not be the main driver of planet formation. Clearly, 
more statistics is needed to obtain robust conclusions. In this 
framework, we report the detection of two planetary companions 
around HIP 11952 as a result of our RV planet survey around 
metal-poor stars. 
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Fig. 1. Comparison of the spectrum of HIP 11952 with stellar 
spectra of the Indo-US library. The dashed lines highlight some 
iron lines (Fe n, Ti n) which are sensitive to the luminosity class 
of the star. The solid lines in the upper part of the figure indicate 
the position of some iron and calcium lines which do not vary 
with the luminosity class. 



The paper is organized as follows: Observations and data 
reduction are presented in section [2] The stellar parameters of 
HIP 11952 are shown in section [3T| together with the most rel- 
evant information available for this star in the literature. In sec- 
tion [X2]we describe the RV and photometric analysis. The de- 
tection of the planetary companion is addressed in section |4] 
Discussion and conclusions are given in sections [5] and [6j re- 
spectively. 



2. Observations and data reduction 

We observed HIP 11952 from August 2009 until January 2011 
with FEROS at the 2.2 m Max-Planck Gesellschaft/European 
Southern Observatory (MPG/ESO) telescope, located at ESO- 
La Silla observatory, during the MPG guarantee time. FEROS 
has a resolution of R = 48 000 and a wavelength coverage of 
370-920 nm ( |Kaufer & Pasquini 1998| . The long-term RV pre- 
cision of FEROS was measured across a period of 6 years, from 
December 2003 until January 2010, using the standard star t 
Ceti (HD 10700). We obtained an RV precision of better than 1 1 
ms _1 . 

The data reduction was performed with a software package 
in the ESO-MIDAS environment available online at the tele- 
scope. The procedure of the RV computation from the FEROS 
spectra is described in Seti awan et al. (2003)| and based on a 
cross-correlation technique. The spectra of HIP 11952 were 



cross-correlated with a template of an F dwarf star which best 
matches the spectral type of HIP 1 1952. 



3. Analysis 

In this section we present the analysis of the stellar parameters 
and RVs derived from the FEROS spectra. 



3.1. The star HIP 11952 



HIP 11952 (HD 16031; LP 710-89) was previously clas 
sified as an F0 dwarf 
Kharchenko & Roeser 2009) 



Wright et al. 2003 



giant G8 (e.g. 



star (e.g., 
and as a 

Sanchez-Blazquez et al. 2006 ICenarro et al. 2007b 

Our spectral classification was carried out by comparing the 
FEROS spectrum of HIP 11952 with spectra from the Indo- 
US library (IValdes et al. 2004i l of metal-poor stars with differ- 
ent luminosity classes and spectral types. The FEROS spectrum 
was convolved to the resolution of the Indo-US spectral library 
(1A). Following the spectral classification criteria presented by 
|Gray & Corbally (2009)| we used some lines (Fe n, Ti n) which 
are sensitive to the luminosity class of the star as well as other 
iron (Fe i) and calcium (Ca i) lines which do not vary with the lu- 
minosity class (Fig.Q]). From this comparison we concluded that 
HIP 11952 is an F2V star. However, this result has to be con- 
firmed by an independent spectroscopic analysis of the stellar 
spectra, as we present below. 

HIP 11952 is at a distance of 115.3 pc as derived from 
the parallax measurements given in the Hipparcos catalogue 
(Perryman 1997). Fundamental parameters of this star were de- 
termined using our high-resolution FEROS spectra. In particu- 
lar, stellar abundances, effective temperature and surface gravity 
were computed using the synthetic spectra from the ID ATLAS 
models ( IKurucz 19931 IKurucz 20051) and the fit of the Ha line 
to the C05BOLD 3D model atmosphere dCaffauet al. 201 11 1. 
The atomic data of the iron lines are from the Large Program 
"First Stars" lead by R. Cayrel, optimized for metal-poor stars 
(see lSivarani et al."2 004). The first attempt of abundance analy- 
sis was based on ID ATLAS model atmospheres computed using 
the version 9 of the co de ATLAS (IKurucz 19931 IKurucz 20531) 
running under Linux dSbordone et al. 20041 ISbordone2 005). 
We derived a temperature of 5960 K by fitting the Ha wings with 
a grid of synthetic spectra computed from ATLAS models, and a 
temperature of 6120 K when we used a grid of synthetic spectra 
based on 3D models. 

By imposing an agreement between the iron abundance de- 
rived from the lines of Fei and the lines of Fen, we derived 
surface gravities log g of 3.8 and 4.0 for the two cases of 
T_eff = 5960 K and 6120K, respectively. The uncertainty on 
the effective temperature derived from the fit of the Ha is 150 K, 
while the error in the surface gravity is 0.3. A microturbulence 
of 1.4 kms -1 was derived by minimizing the slope of the abun- 
dance versus equivalent width (EW) relation. The resulting iron 
abundance derived is [Fe/H]= -1.95 + 0.09 for T eff = 5960 K 
and [Fe/H]= -1.85 ± 0.09 for T eff = 6120 K, respectively. 

We thus obtained two possible parameter sets for this 
star: (T eff , logg, [Fe/H])=(5960 K, 3.8,-1.95) (ID ATLAS) 
and (6120 K, 4.0, o -1.85) (3D MODELS). We measured an 
EW=3.07+0.03 mA of the Li feature at 670.7 nm, which implies 
an abundance A(Li)=2.2 if we fix T e ff = 6120 K, and A(Li)=2.1 
in the case T e ff = 5960 K. 

Feltzing et al. (2001) derived a stellar age of 12.8 + 2.6 Gyr 



by comparing the Stromgren photometry of HIP 11952 with 



2 



J. Setiawan et al.: Planetary companions around the metal-poor star HIP 1 1952 



the evolutionary tracks computed for the Stromgren metallic- 
ity ([m/H]= -1.54) of the star. We checked their result using 
the photometry provided by Hipparcos for HIP 11952 in the 
Bessell filters system (Bessell 2000) with the isochrones by 
Marigo et al. (2008)| and |Girardi et al. (2010)1 calculated for the 



Stromgren metallicity of the star. In the assumption that the dust 
extinction along the line of sight to HIP 11952 is not large, the 
Hipparcos photometry indicates, within its uncertainty, an age 
older than 10 Gyr for this star. 
We 



Marigo et al. (2008) and 



used the isochrones by 
|Girardi et al. (2010)| computed for [m/H]= -1.54 and ages in the 
range 10-13 Gyr in order to consistenly derive the mass and ra- 
dius of HIP 11952. By comparing its photometry with the se- 
lected isochrones, we constrained its mass between 0.79 M and 
0.88 M . From the relation between surface gravity, stellar mass 
and stellar radius: 



(1) 



we derived a stellar radius of 1 .6+0. 1 R , where the error is ob- 
tained from the errors propagation. The derived stellar mass and 
radius adopted in this work are: 0.83^4 M and 1.6+0.1 R . 

Fundamental stellar parameters have also been es- 
timated in previous studies (e.g., Ma sana et al. 2 006 ; 
Char bonnel et al. 20 05 ). The corresponding values are 
T e ff = 6367 K, a stellar radius R* = 1 .0 R Q and surface gravity 




4.1. |Feltzing et al. (200 1)1 reported m = 0.785 M . If 
we compile all literature values for the surface gravity given 
Cayrel de Strobel et al. (200T)| and use R, = 1.0 

± 0.23 M 



, we 

Within the errors, 



obtain a mean value m = 0.55 
these values are consistent with those we derived, although the 
stellar radius is smaller than the one we adopted. Finally, we 
compared the derived stellar parameters with those given by 



Casagrande et al. (2010) These authors measured log g =4.17, 
[Fe/H]=-1.74 and T eff = 6186 K. Within the uncertainties, 
these values are in good agreement with our determination. 

Based on the surface gravity and radius derived from our 
analysis and the available literature data, HIP 11952 is more 
likely a star already evolved off the main-sequence, roughly sit- 
ting at the base of the subgiant branch. 

Astrometric and photometric data of HIP 11952 can 
be found in public catalogues (e.g., Hipparcos). The as- 
trometric variability is less than 5 mas, allowing the 
conclusion that an unseen stellar companion in the sys- 
tem can be ruled out. Nevertheless, older RV measure- 
ments from 



Eggen & Sandage (1959) 



combined with 
HIP 11952 may 
However, our 



Carney & Latham (1987) suggested that 
be a spectroscopic binary SB1 dFouts 198Tb 
measurements can neither reject nor confirm this claim yet. The 
updated parameters of HIP 1 1952 are given in TableQ] 

HIP 11952 is listed as a member of the metal-poor stellar 
stream detected by |Arifyanto & Fuchs (2006)| a group of pu- 
tative thick disk stars moving on similar orbits in the Galactic 
potential and currently lagging the Local Standard of Rest by 
Viag ~ 80 km s . As such, HIP 11952 might stem from one of 
the Milky Way's former satellite galaxies that once got tidally 
disrupted. Alternatively, Mi nchev et al. (2009)| showed that the 
stream might consist of in situ disk stars being scattered towards 
common orbits through a merger-induced perturbation of the old 
disk ~ 1 .9 Gyr ago. It is also possible that there is a connection 
to the Arcturus stream at Vi ag w 100 which itself could have a 
tidal or a dynamical origin (see discussion in Kle ment (2010)| for 
more details). 




5200 5400 

time (JD-2450000) 




Period (d) 

Fig. 2. RV observations of HIP 1 1952, taken with FEROS from 
August 2009 until January 2011 (upper panel). We calculated a 
single Keplerian fit to the data (solid line) for the outer compo- 
nent. A multi-component Keplerian fit is presented in Fig|9] see 
Table. [3] The lower panel shows a Generalized Lomb-Scargle 
(GLS) periodogram of the whole RV data set. The highest peak, 
marked with Pi corresponds to a period of 290 + 16 d with a 
False Alarm Probability of ~ 7 x 10~ 6 . 

Table 1. Stellar parameters of HIP 1 1952 



Parameter 


value 


unit 


Spectral type* F2V - IV 




(a) 

my 


9.88 ± 0.02 


mag 


Parallax'"' 


8.67 ± 1.81 


mas 


T* 

1 ell- 


5960 ± 150 


K 


6120 ± 150 


K 


log g* 


3.8 + 0.3 


cm 2 /g 
cm 2 /g 




4.0 + 0.3 


[Fe/H]* 


-1.95 + 0.09 


dex 




-1.85 + 0.09 


dex 


r: 


1.6 + 0.1 


Ro 


Mass* 


0.833E 


M G 


Age<*» 


12.8 + 2.6 


Gyr 


Vrot sin i* 


5.2+1.0 


km s -1 


P ro ,/ sin i 


15.7 + 2.5 


d 



* this work 

" Hipparcos catalogue ( Perryman 1997[ > 
A |Feltzingetal. (2001)| 



3.2. Radial velocity 



The RV variation of HIP 1 1952 is shown in Fig.|2](upper panel). 
During the observation campaigns we obtained 77 RV mea- 
surements (Table. 13. We applied the Generalized Lomb-Scargle 
(GLS) periodogram dZechmeister & Ktirster 20091 ) to the RV 
data in order to search for periodicities. 

We found several signals in the periodogram, as shown in 
the lower panel of Fig. [2] The four highest signals are marked 
with Pi, P2, P3 and P4. The highest peak Pi has a False Alarm 
Probability (FAP) of 7.2 x 10~ 6 and corresponds to a period of 
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Fig. 3. The residual RV after removing the 290 d periodicity. The 
periodogram in the upper panel shows two peaks at 6.95 d and 
1.16 d. The 1.16 d is identified as a harmonic of 6.95 d (1/1.16 + 
1/6.95 = 1.0). The lower panel shows the phase folded RVs with 
a period of 6.95 d. The solid line shows a Keplerian fit for the 
residual RV variation. 



290 + 16 d. We computed a Keplerian fit to the data, shown as 
the solid line in the upper panel of Fig. [2] The parameters of this 
fit are given in Table [3] 

After removing the 290 d signal, the peaks P2 and P3 disap- 
pear. However, the peak P4 remains, which means that this sig- 
nal is not an alias of P\, The signal P4 corresponds to a period 
of 6.95 d with FAP= 8 x 10 4 . In the residual RV periodogram, 
we also observed a signal at 1.16 d (Fig. [3]upper panel), which 
is obviously a harmonic of the 6.95 d period (1/6.95 + 1/1.16 
= 1.0). In the lower panel of Fig. [3] we have phase folded the 
RVs with P = 6.95 d and show a Keplerian fit with the parame- 
ters listed in Table[3] After excluding stellar activity as the cause 
for the RV variations in the next sections, we are going to inter- 
pret the signals P\ and P4 as orbital periods of unseen low-mass 
companions (section|4]below). 



Pj sin/ = 2nR*lvsmi, the maximum stellar rotation period is 
found to be 15.7+2.5 d. 



3.3.1. Line profile asymmetry 

The line profile asymmetry can be quantified by the bisector ve- 
locity span (BVS). A definition of the BVS is given, e.g., in 
Hatzes (1996) We measured the BVS of the stellar spectra and 
searched for any periodicity that might be related to the RV vari- 
ation. In the GLS periodogram of BVS we found, however, no 
significant peak. Thus, we cannot use the bisector to determine 
the stellar rotation period. We then searched for a correlation 
between BVS and RV to find out whether the observed RV vari- 
ation is due to rotational modulation. We computed a correlation 
coefficient c — 0.1 between the RV and BVS (Fig. |4j». This value 
indicates that the RV variation is not correlated with the BVS. 
However, it does not give any hint about the stellar rotation pe- 
riod. 



-0.5 




RV (km s -1 ) 

Fig. 4. The figure shows the measurements of BVS, ploted 
against residual RVs. The plot shows no correlation between 
BVS and RV. We computed a correlation coefficient c — 0.1. 



3.3. Stellar rotation 

A systematic investigation of the stellar activity is mandatory 
to avoid wrong interpretations of the observed RV variations. 
There are several possibilities to probe the stellar activity. The 
line profile asymmetry (bisector) and the Ca n lines are known 
as reliable stellar activity indicators. These activity indicators, 
if they show periodic variations, can be used to determine the 
stellar rotation period. Besides the spectroscopic methods, pho- 
tometric observations are also commonly used to find the stellar 
rotation period. 

HIP 11952 itself is not a star with high stellar activity. This 
is inferred from the absence of emission cores in Ca 11 K (/13934) 
and H (/13967). Furthermore, no Ha emission line was observed 
in the spectra. The projected rotational velocity is also rela- 
tively low (vsin; = 5.2 km s" 1 ) compared to other F-type dwarf 
stars ( de Medeiros & Mayor 1999} . Finally, using the relation 



3.3.2. Can analysis 

As mentioned before we found no emission cores in the Ca n H 
& K lines (/1/13934,3968) which are, in general, excellent stel- 
lar activity indicators to probe the rotational modulation of the 
star. Nevertheless, we investigated to possibility to use Ca 11 K 
lines to find any indication of stellar activity. We calculated the 
stellar activity index, known as S -index, similar to the method 
described in |Vaughan et al. ( 1978)] and |Santos et al. (200"0)1 We 
found a significant periodicity of 36 d in the periodogram, but 
the error bar of the individual measurement is large. This value 
is also close to the typical observational window of about one 
month. Thus, also by considering Pj sin i value, we do not adopt 
this as the stellar rotation period. 

We exploited the capability of FEROS to investigate the Can 
lines /1/18498,8662, following the technique presented e.g. by 
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Table 2. RV variation of HIP 1 1952 



JD RV error 

-2450000 (m/s) (m/s) 
5052.8320 24262.99 46.17 
5053.9449 24398.12 26.57 
5054.9025 24388.69 26.20 
5056.7834 24447.32 31.11 
5057.9442 24257.13 20.34 
5060.8739 24280.16 21.22 
5063.8340 24418.54 41.10 
5064.8816 24338.95 28.32 
5068.8975 24302.92 34.02 
5071.9278 24242.14 22.27 
5072.9043 24306.86 24.57 
5073.7695 24300.97 24.45 
5074.9106 24317.19 28.81 
5075.8276 24375.46 30.57 

5107.7807 24126.26 44.43 
5108.8472 24079.34 39.29 
5109.7623 24143.82 30.85 
5109.7800 24221.74 27.24 
5110.7724 24281.43 29.66 
5112.7689 24325.34 28.45 

5112.7808 24281.70 29.68 
5113.7627 24226.34 32.90 
5164.7436 24276.41 29.50 
5166.7405 24373.37 34.20 
5167.6315 24342.35 32.31 
5168.6411 24288.16 21.91 
5170.7375 24230.07 34.85 
5172.6838 24297.89 41.50 
5174.6691 24319.12 26.32 
5181.5403 24428.12 32.69 
5198.5285 24328.64 35.72 
5199.6074 24340.87 28.77 
5200.6219 24390.29 30.19 
5201.5849 24400.51 21.81 
5202.6165 24441.51 25.55 
5203.5918 24288.39 30.66 
5204.6389 24263.70 31.88 
5254.5427 24299.59 41.61 
5255.5272 24402.04 40.17 



JD RV error 

-2450000 (m/s) (m/s) 

5256.5331 24342.01 35.92 

5258.5279 24461.25 33.87 

5262.5458 24312.08 39.98 

5268.5034 24367.11 33.39 

5349.9353 24151.48 35.26 

5351.9288 24351.42 26.78 

5357.9251 24358.36 50.76 

5358.9292 24410.18 36.79 

5396.8522 24319.13 45.60 

5399.8478 24342.35 40.48 

5400.9186 24334.94 34.61 

5401.8807 24423.26 29.23 

5402.8758 24473.63 38.74 

5403.8058 24401.82 37.24 

5404.8730 24261.93 41.39 

5405.8598 24302.11 35.31 

5407.8695 24314.65 34.22 

5430.8126 24397.33 32.08 

5432.7291 24217.71 40.97 

5434.8370 24277.26 30.19 

5436.8532 24391.30 29.11 

5437.8192 24431.76 34.17 

5491.5926 24425.56 46.26 

5492.7290 24320.54 42.54 

5493.7936 24351.56 32.04 

5494.8416 24371.35 37.64 

5497.8011 24298.56 30.49 

5521.6567 24476.62 22.33 

5523.6496 24471.39 25.01 

5524.5791 24387.17 17.30 

5525.6558 24390.15 17.25 

5526.6492 24439.61 15.05 

5527.6314 24449.32 25.68 

5583.6321 24489.01 14.01 

5584.5958 24630.11 21.68 

5585.6161 24437.91 36.76 

5587.6035 24290.41 50.00 

5588.6275 24391.65 66.95 



|Larson (19 93), where Ca n /18662 was used to determine the stel- 
lar rotation period. [Setiawan et al. (2010)] for example, used the 
equivalent width (EW) variations of Ca n /18498 as a stellar ac- 
tivity indicator to estimate a stellar rotation period that agrees 
with the bisector analysis. 

The EW variation of Ca ,18498 of HIP 1 1952 indeed shows 
a periodic variation with P - 1 .76 d with a FAP of few percent. 
Thus, it is only marginaly significant. Following |Larson (1993)| 
we then measured the EW variation of Ca n /18662. Interestingly, 
we found a significant periodicity in the EW variation of Can 
T8662. The signal corresponds to a peak at a period of P — 4.82 
d with FAP= 8 x 10 3 . Fig.|6]shows the EW variation and GLS 
periodogram of the Ca n /18662 line. Assuming that this feature 
is related to the stellar magnetic activity caused by starspots, the 
period is most-likely linked to the stellar rotation. 

3.3.3. Photometric data 

Photometric V band observations of HIP 11952 are available 
in the Hipparcos catalogue. Unfortunately, the data set is very 
sparse, with 72 photometric measurements over a time span 



of 893 d. The minimum time-interval between data points is 
a;0.014d, and the typical photometric errors are ^0.02 mag. 
Due to long-term gaps of several days in the data set, period- 
icities of a few days cannot be reliably detected. The sampling 
allows for the detection of very short-term (few hours) as well as 
long-term (20 d) variations. 

We searched for periodicity in the photometric data us- 
ing a combination of two periodogram analysis techniques: 
the Scargle periodogram ( Scargle 19 82) and the CLEAN 
algorithm (Rober ts"1987l ). The combination of these two 
techniques provides a reliable period detection as outlined 
in several rotational period and variability studies (e.g., 



Rodriguez-Ledesma et al. 2009 1. 



Based on the Scargle periodogram, we detected significant 
signals at P l = 0.072, P 2 = 0.33 and at P 3 =2 d. When 
the CLEAN algorithm is applied, the 2 d signal is removed 
and therefore, we concluded that it is probably a false peak 
or alias due to the clumpy data sampling. Both the 0.072 and 
0.33 d periods in the power spectrum remains (Fig. [7]), with 
FAP, based on the Scargle periodogram, of 0.5 % and 1 %, 
respectively. All other peaks in Fig. 6 have larger FAPs. We 
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Fig. 5. Upper panel: The measurements of Ca n /13934 S -index. 
Lower panel: The GLS periodogram of Ca n /13934 shows a clear 
peak corresponding to a period of P — 36 d. However, the error 
bar of the individual measurement is large and the period is close 
to the observational time window of about 1 month. 



respectively. Fig.[8]shows the phase folded light curves with the 
periods of 0.072 and 0.33 d. 

Due to the quality of the data set, however, it is difficult to 
ensure the significance of these periods. Our analysis suggests 
some evidence for short-term photometric variations, which 
might be interpreted as possible pulsation modes in this F-type 
star. 



< 

LU 
_1 

o 




1 10 
Period [day] 



Fig. 7. The photometric observations of HIP 11952 show two 
peaks at P = 0.3 and P = 0.072 d. 
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Fig. 6. Upper panel: The measurements of Ca n /18662 EW vari- 
ation. Lower panel: The GLS periodogram of Ca n /18662 shows 
a clear peak corresponding to a period of P — 4.82 d. The dashed 
line shows the period of the RV variations. 



have also computed the statistical F-test and a derived FAP 
from it (Scholz ~et al. 2004l and Rodriguez-Ledesma et al. 2009). 
The FAPf, fs/ represents the probability that the period found is 
caused by variations in the photometric noise, and therefore it 
is independent of the periodogram analysis. The FAP Ftest de- 
rived for the 0.072 and 0.33 d detected peaks are 5 % and 11%, 




Phase 
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Fig. 8. Phase folded light curves from the Hipparcos photometric 
data. 



Based on the analysis of the stellar activity indicators, we as- 
sume a stellar rotation period of 4.8 d, as derived from the Can 
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/18662. Additionally, there might be evidence for stellar pulsa- 
tions with intra-day periodicities 

4. Planetary companion 

Since the long-period and short-period RV variations have dif- 
ferent characteristics from those of the stellar activity indicators, 
we concluded that they are most likely caused by the presence 
of unseen companions. 

We computed the orbital solution by using a two- 
components Keplerian model. In Fig. [9] we show the calcu- 
lated orbital fit and residual velocities. The orbital paramaters 
are given in Table [3] With a derived primary mass of 0.83 
M , we calculated the minimum masses of the companions 
m.2 sin i = 0.78 Mj up for the inner and 1112 sin i = 2.93 Mj up for 
the outer companion. The orbital semi-major axes are 0.07 AU 
and 0.81 AU, respectively. The planetary orbits have moderate 
eccentricities of 0.35 and 0.27, which seem to be not unusual, 
based on the statistics of the eccentricity of exoplanets (see e.g., 
www.exoplanet.eu). 




time (JD-3450DD0) 



Fig. 9. Orbital solution of a two-planets Keplerian fit (upper 
panel). The orbital parameters are listed in Tablef3] The resid- 
ual velocities are shown in the lower panel. 



Table 3. Orbital parameters for HIP 1 1952 b and c 



Parameter 


Unit 


HIP 11952b 


HIP 11952c 


P 


d 


290.0 ± 16.2 


6.95 ± 0.01 


T 


JD 


5402.0+ 1.3 


5029.2 ± 0.04 




-2450000 






e 




0.27 ±0.10 


0.35 ± 0.24 


w. 


deg 


59.3 ± 2.5 


61.2 ±6.6 


Ki 


ms" 1 


105.2 ± 14.7 


100.3 ± 19.4 


mi sin i 


M Jup 


2.93 ± 0.42 


0.78 ±0.16 


a 


AU 


0.81 ±0.02 


0.07 ± 0.01 


V 


kms" 1 


24.365 ± 0.01 




a(0 - C) 


ms" 1 


70.22 




reduced^- 2 




1.3 





When calculating the orbital solutions, we obtained a rela- 
tively large cr(0 - C) value. A possible explanation to this is the 



presence of another unseen companion or RV jitter due to the in- 
teraction between the two companions. Because of the absence 
of the emission cores in Ca 11 H & K, the large cr(0 - C) is most 
likely not due to the intrinsic stellar variability. 

We investigated the residual velocities and found a signifi- 
cant signal at ~40 days. However, the amplitude of the residual 
RV variations is in the order of the error bars. Moreover, the 
window function shows also a peak close to this period. To de- 
tect further low-amplitude RV variations, more intensive high- 
precise RV measurements are needed. 

From the knowledge of P ml / sin i = 15.7 d (TableQ]) and the 
rotation period P = 4.82 d, we derived an inclination angle of 
the stellar rotation of 18°. Assuming that the orbital inclination 
of the companion does not differ much from the stellar rotation 
inclination angle, we estimated true companion masses of 2.5 
and 9.5 M Jup . 

5. Discussion 

A fundamental parameter of HIP 11952 is the stellar metallic - 
ity. The metallicity issue here is particularly interesting, since 
according to the theory of planet formation via core-accretion 
processes, planetary companions around such metal poor stars 
like HIP 1 1952 are not expected. The majority of the planet host 
main-sequence stars are metal rich (Fig. [TO] upper panel). 
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Fig. 10. The metallicity distribution of stars hosting planets. 
Solid lines in the top and bottom panels represent dwarfs 
and giants, respectively (data from the exoplanet encyclopedia 
www.exoplanet.eu). Dash lines show data from the exoplanet 
orbit database (Wri ght et al. 201 ty which does not discriminate 
between dwarf and giant stars. 

For dwarf stars, the metallicity generally reflects the ini- 
tial metallicity at star formation. Assuming that HIP 11952 is 
a dwarf or a turn off star, its initial low metallicity makes 
HIP 11952 unusual among the planetary systems discovered so 
far (Fig. ED. 
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The both planetary companions around HIP 11952 belong 
to only few planets that have been discovered in low metallic- 
ity systems ([Fe/H]<-1.5), comparable to HIP 13044 reported 
by |Setiawan et al. (2010)[ Note that, so far also only few plan- 
ets have been detected with host star's metallicities -l<[Fe/H]<- 
0.5. This group includes -60% main-sequence and ~40% giants. 
From the current statistics, about 50% of the giant planet-host 
stars have sub-solar metallicities. 

Whether the metallicity of a giant reflects its initial metal- 
licity is still under debate. The convective envelopes in giant 
star are much more extended than in main-sequence stars (e.g. 
Pasquini et al. 2007) and they might alter their surface chemical 
abundances. 

The presence and formation of planets around metal-poor 
stars, in particular those with metallicities [Fe/H]< -1.5 are 
still poorly understood. According to the core accretion the- 
ory (e.g., ISafronov 19691 IPollack et al. 1996t high metallicity 
is required for the formation of planets. Alternatively, planets 
around metal-poor stars could form by gravitational disk in- 
stability processes (e.g., IBoss 1998b or other mechanisms (see 
e.g., Nayakshin 2011). According to Nayakshin (2011) planets 
around such metal-poor stars have no solid cores and may form 
as a result of the second collapse (H2 disassociation) of their 
embryos and radial migrations. Further discoveries of planets 
around metal-poor stars can provide more constraints on cur- 
rently different planet formation theories. 

The planets around HIP 11952 have orbits with moderate 
eccentricities. To examine the dynamical stability of the sys- 
tem, numerical simulations are necessary to find stable planetary 
configurations, see e.g., Barnes & Quinn (2004) Following their 
calculations, systems with two or more planets in large separated 
orbits are fully stable. The orbits of HIP 1 1952 b and c are in a 
42: 1 ratio and thus far beyond the 10: 1 resonance. Therefore, the 
system HIP 1 1952 is most likely fully stable. 

Finally, HIP 1 1952 and its planets are among the oldest plan- 
etary systems known. HIP 1 1952 is also older than typical stellar 
ages in the Galactic thick disk. A possible connection to a metal- 
poor stellar stream reported by Arify anto & Fuchs (2006)| is an 
interesting aspect since it suggests that HIP 11952 might actu- 
ally belong to a part of the thick disk that has been accreted from 
a disrupted former satellite galaxy, similar to HIP 13044. The 



age estimation of 12.8 ± 2.6 Gyr given by Feltzing et al. (2001) 



is 1 Gyr older than that of HD 114752 which has an estimated 
age of 11.8 + 3.9 Gyr. HIP 1 1952's age is close to the one of 
HE 15 23-0901, which is the oldest star (13.2 + 2.7 Gyr) known 
so far (IFrebel et al. 2 007). The old stellar age is supported by the 
very low metallicity of HIP 1 1952 which is typical of Population 
II stars. 

6. Conclusions 

We observed RV variations of HIP 1 1952. The spectroscopic and 
photometric analysis of the star show that the periodic RV vari- 
ations are not caused by the intrinsic stellar variability. Based 
on our analysis, we detected two planetary companions around 
the metal poor star HIP 1 1952 with orbital periods of 6.95 d and 
290 d. We found evidence for intra-day stellar pulsations and ob- 
served a stellar rotation of 4.82 d. We computed the companion's 
minimum mass of mi sin / = 0.78 Mj up for the inner planet and 
m.2 sin i = 2.93 Mj up for the outer one. Additional high-precise 
RV measurements are necessary to improve the orbital solution 
and put more constraints on the eccentricities. Further RV obser- 
vations might also reveal the presence of other low-mass com- 
panions in the system. From the metal abundance analysis that 



we carried out, we obtained an average [Fe/H]= -1.90 + 0.06 
from Fe i and Fe n, respectively, which makes HIP 1 1952 b and 
c the first planets discovered around a dwarf or subgiant star with 
[Fe/H]<-1.5. This discovery is also remarkable since the plane- 
tary system most likely belong to the first generation of planetary 
systems in the Milky Way. 
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